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ABSTRACT 



INTRODUCTION 



The usefulness of studying elliptical galaxies comes both 
from the importance for cosmology of knowing their 
past luminosities and colours, and for testing stellar evo- 
lution models and theories on supernova progenitors. 

Many features shown by elliptical galaxies required, 
in the past years, theoretical explanations. For example, 
in order to reproduce the well known relation between 
mass and metallicity Mathews and Baker (1971) and 
Larson (1974) proposed the existence of galactic winds 
driven by supernova explosions (note that with such a 
way they explained also the apparent lack of gas in these 
systems) . 

Successively, a great deal of effort has been de- 
voted in studying detailed chemical evolution of ellipti- 
cal galaxies (Arimoto and Yoshii, 1986, 1987; Matteucci 
and Tornambe, 1987; Brocato E., Matteucci F., Mazz- 
itelli I., Tornambe A., 1990; Angeletti and Giannone, 
1990; Matteucci, 1992; Ferrini and Poggianti, 1993; Mat- 



Models of chemical evolution of elliptical galaxies taking into account dif- 
ferent escape velocities at different galactocentric radii are presented. As a 
consequence of this, the chemical evolution develops differently in different 
galactic regions; in particular, we find that the galactic wind, powered by 
supernovae (of type II and I) starts, under suitable conditions, in the outer 
regions and successively develops in the central ones. The rate of star forma- 
tion (SFR) is assumed to stop after the onset of the galactic wind in each 
region. The main result found in the present work is that this mechanism is 
able to reproduce metallicity gradients, namely the gradients in the Mg^ in- 
dex, in good agreement with observational data. We also find that in order to 
honor the constant [Mg/Fe] ratio with galactocentric distance, as inferred from 
metallicity indices, a variable initial mass function as a function of galactocen- 
tric distance is required. This is only a suggestion since trends on abundances 
inferred just from metallicity indices are still uncertain. 

teucci, 1994; Bressan, Chiosi and Fagotto, 1994; Mat- 
teucci and Gibson, 1995; Gibson, 1997; Gibson and Mat- 
teucci 1997). For example, on the basis of supernova 
driven galactic wind models, Matteucci and Tornambe 
(1987) first analyzed the evolution of Mg and Fe in these 
galaxies and predicted that [Mg/Fe] should decrease as 
a function of the galactic mass, at variance with what 
has been more recently suggested by population syn- 
thesis studies (Worthey, Faber and Gonzalez, 1992) in- 
dicating that [Mg/Fe] increases with galactic mass in 
the nuclei of giant ellipticals. Matteucci (1994) pointed 
out that, in order to explain the origin of this relation 
one has to assume either that the efficency of star for- 
mation is an increasing function of mass, or that the 
initial mass function (IMF) should favor more massive 
stars in more massive galaxies. In the models of Larson 
(1974), Arimoto and Yoshii (1987), and Matteucci and 
Tornambe (1987), in fact, the star formation efficiency 
was either assumed to be constant or to decrease with 
galactic mass. This last hypothesis was based upon the 
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fact that star formation is mainly determined by cloud- 
cloud collisions. In this frame-work the time scale for 
star formation (namely the time required to consume 
all the gas in a galaxy by star formation and corre- 
sponding to the inverse of the star formation efficiency) 
was assumed to be proportional to the cloud-cloud col- 
lision time or to the dynamical free-fall time, with the 
consequence of decreasing with the galactic mass. 

Concerning the other explanation relative to the 
IMF, the observed variation of the M/Lb ratio in ellipti- 
cal galaxies as a function of their blue luminosity (Faber 
and Jackson, 1976; Michard, 1983; Kormendy and Djor- 
govski, 1989; van den Marel, 1991; Bender, Burstein and 
Faber, 1992) seems to argue in favor of a possible vari- 
ation of the IMF with galactic mass, although other 
possible explanations cannot be excluded. Renzini and 
Ciotti (1993), for example, proposed a variable amount 
and/or concentration of dark matter as a function of lu- 
minous mass; Zepf and Silk (1996) analyzed the hypote- 
sis that star formation includes a phase in which star 
formation is restricted to massive stars, with the bias 
toward high mass stars increasing with galactic mass. 

The evidence for a mass-metallicity relation in el- 
lipticals comes not only from changes in colors with 
galactic luminosity (Faber 1972), but also from the lin- 
ear relationship existing between the Mg^ index as de- 
fined by Faber S.M., Friel E.D., Burnstein D., Gaskell 
CM. (1985), and the internal velocity dispersion a 
(Dressier A., Lynden-Bell D., Burnstein D., Davies R.L., 
Faber S.M., Terlevich R.J. and Wegner G. 1987, Ben- 
der, Burstein and Faber, 1993). It is generally assumed, 
in fact, that these color and metallicity index changes 
indicate variations in metallicity. 

Ziegler (1997) found that determining the Mg b — a 
(Mgb is another metallicity indicator, always referring 
to Mg, as defined in Faber et al. 1985) relations at var- 
ious redshifts is a powerful and robust method to mea- 
sure the evolution of elliptical galaxies. In his work he 
points out two conclusions: first, the stellar population 
of elliptical cluster galaxies evolves passively between 
z = 0.4 and today, and second, the epoch of formation 
of the stars which form the present day ellipticals oc- 
curred at high redshifts( Z > 2 and the most luminous 
objects might have formed even at Z > 4). 

Another important aspect of the Mgb — a relation 
is its small scatter. This is usually interpreted as an in- 
dication that the process of galaxy formation was short 
(see Renzini 1995) 

Metallicity gradients are also well known charac- 
teristics of elliptical galaxies (Carollo Danziger and Bu- 
son, 1993; Davies, Sadler and Peletier, 1993) The best 



evidence comes from the spectroscopic measurements 
of absorption lines in ellipticals (Faber 1977; Burnstein 
D., Faber S. M., Gaskell C. M., Krumn N, 1984; Faber 
et al. 1985; Peletier 1989). Typically, the Mg b and 
Fe absorption-line strengths are stronger in high-mass 
galaxies. In individual galaxies they generally decrease 
outward from the center. The colors behave similarly, 
being redder in high-mass galaxies and changing from 
red to blue outwards. Therefore, both the color and 
line strength changes indicate consistently that these 
gradients are very likely due to variations in metallic- 
ity. However, the mechanism by which these metallicity 
variations could have arisen is still unclear. The models 
proposed are based on the processes developed during 
galaxy formation. Larson (1974) and Carlberg (1984) 
showed that elliptical galaxies which form by dissipative 
collapse of a primordial gas gas cloud, should have an in- 
creasing metallicity with decreasing radius. The idea is 
that the stars begin to form everywhere during the col- 
lapse and remain in orbit with a little inward motion, 
whereas the gas falls towards the center and is enriched 
by the evolving stars. In this way, the stars which form 
at the galaxy center are more metal enriched than those 
formed in the external regions. These models produce 
galaxies with very steep metallicity gradients: a factor 
of 10 reduction in metallicity for a factor of 10 increase 
in radius. Carlberg's models include a pressure term as- 
sociated with the energy input from supernovae during 
the collapse. This decreases the size of abundance gra- 
dients in the final galaxies to a factor of 3 over each 
decade in radius. 

Franx and Illingworth (1990) found that, for their 
sample of 17 galaxies, the local (B-R) (U-R) colors are 
functions of the local escape velocity v eac for all of 
the galaxies. This result was confirmed later by Davies 
Sadler and Peletier (1993) and Carollo and Danziger 
(1994) using a more direct indicator of metallicity 
(Mg?) than colors and a more appropriate calculation of 
v esc for each galaxy. This relation suggests that metal- 
licity gradients could have arisen because of the different 
times of occurence of galactic winds in different galactic 
regions. In fact the galactic wind starts when the energy 
injected in the interstellar medium (ISM) by supernovae 
(of both type I and II) becomes equal to the binding 
energy of gas (Larson 1974). Moreover, it is reasonable 
to think that, when the galactic wind starts, the SFR 
breaks down. Hence, in the regions where low (i.e. 

where the local potential is shallow), the galactic wind 
will develops earlier and the gas is less processed than 
in the regions where v esc is higher, and the potential 
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well in the external regions is lower than in the internal 
ones. 

The aim of this paper is to discuss the forma- 
tion of abundance gradients in elliptical galaxies. The 
model suggested here to explain metallicity gradients 
does not follow in detail the processes occurred during 
galaxy formation such as collapse or merging, and it 
starts with the total galactic mass already present at 
the beginning. This is a good assumption for elliptical 
galaxies since they must have formed in a reasonably 
short timescale (Matteucci 1994, Renzini 1995, Ziegler 
1997), as already discussed before. Recently, Kauffman 
and Chariot (1997) have discussed the formation of el- 
lipticals by merging of spirals in the context of a hierar- 
chical clustering scenario. This mechanism would imply 
a much longer process than that adopted here and it 
certainly can not explain the high [Mg/Fe] ratios found 
in the nuclei of ellipticals. Moreover, recently Prochaska 
and Wolfe (1997) showed that Kauffman (1996) model 
is inconsistent with the damped Lyman-a data. 

The results we present here are obtained using a 
code of chemical evolution (Matteucci and Gibson 1995, 
hereafter MG) modified in order to take account the 
different evolution of each galactic region. 

The work is organized in four sections. Section 2 
concerns the new theoretical prescriptions with respect 
to MG; the theoretical results are discussed and com- 
pared with observations in section 3, where, models of 
a galaxy with a mass of 10 n M© are discussed for illus- 
trative purposes. Finally, some conclusions and remarks 
are done in section 4. 

2 THEORETICAL PRESCRIPTIONS 

2.1 The chemical evolution model 

In the following, we consider elliptical galaxies to have a 
spherical symmetry. We assume that at the initial time 
(t=0) these systems consist only of primordial gas and 
the star formation starts at the same time in all galactic 
regions. 

The galaxy is partitioned in zones having the shape 
of spherical shell. The basic equations to calculate the 
chemical evolution in each zone are the same as in MG 
and, for this reason, they will not be repeated here. We 
will only remind that the model belongs to the category 
of supernova driven wind models and that follows the 
evolution in space and time of the abundances of several 
chemical elements (H, D, 3 He, 4 He, C, N, O, Ne, Mg, 
Si, Fe and others). The stellar lifetimes are taken into 
account and the rates of supernovae of type II and I 
are computed in detail. In particular, for the rate of 



type la supernovae, which are assumed to originate in 
white dwarfs in binary systems we adopt the original 
formulation described in Matteucci and Greggio (1986) 
and Matteucci and Tornambe (1987). 

2.2 Stellar Birthrate Function 

The stellar birthrate function is generally separated into 
two independent functions: 

B(m,t) = ip(t)4>(m) (1) 
Following MG, ip(t) is given by: 

iP(t) = v Pgas {t) (2) 

where p ga s{t) is the volume gas density in the region of 
interest at the time t. The parameter v represents the 
star formation efficiency and is taken from Arimoto and 
Yoshii (1986) v = 8.6(M 9as (0)/10 12 M Q )-°' 115 (Gyr- 1 ), 
and is assumed to be constant with radius. 

The IMF is generally expressed as a power law and 
here we consider: 

4>(m) oc m~ x (3) 

normalized to unity in the mass range 0.1 < M/Mq < 
100. We have explored three different cases: 

(i) x — 0.95 for all the regions of galaxy; 

(ii) x = 1.35 for all the regions of galaxy; 

(iii) x variable and in particular increasing with ra- 
dius. 

2.3 Stellar Yields 

We adopted the same stellar yields as in MG in order to 
be able to compare our results with theirs. In particular 
we assumed: 

i) for massive stars and supernovae II (M > SMq) 
the yields of Woosley (1987) are adopted, 

ii) for low and intermediate mass stars (0.8 < 
M/Mq < 8) we assume the yields of Renzini and Voli 
(1981). 

iii) For type la supernovae the yields of Thiele- 
mann, Nomoto and Yokoi (1993) for their model W7 
are used. The stellar lifetimes are the same adopted by 
Matteucci and Padovani (1993), namely: 

for M > 6.6M 

r m = 1.2m" 188 + 0.003 Gyr (4) 
and for M < 6.6M 

Tm _ 10 1. 338-^1. 790-0. 2232-(7.764-!o 9 m) _ g ^ 

For the details about calculation of the type la SN 
rate we address the reader to MG95. We will remind 
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here only that the fraction of binary systems is a func- 
tion of the assumed IMF as described in MG95. 



2.4 Galactic Wind 

Galactic winds in ellipticals were originally introduced 
to reproduce the well known mass-metallicity relation 
in these systems and to explain the apparent lack of gas 
(Larson 1974). However, besides this, galactic winds are 
the natural consequence of a phase of intense star forma- 
tion. The existence of a galactic wind, in fact, is strictly 
related to the energetics of the interstellar medium com- 
pared to the potential energy of the gas. The time at 
which a galactic wind occurs is a fundamental param- 
eter for the chemical evolution of a galaxy region. The 
condition which should be satisfied for the gas to be 
ejected from a given region is expressed by the follow- 
ing condition (Larson 1974): 



El hsN (taw) = E l Bgas (tGw) 



(6) 



where the index i indicates the particular region under 
consideration. 

The condition requires that the thermal energy of the 
gas in a given region, resulting from supernova explo- 
sions in that region, must exceed the binding energy 
of the gas. For both quantities in equation (6) we used 
new theoretical prescriptions with respect to previous 
works. These prescriptions will be described in the next 
paragraphs. 

2.5 The potential energy of the gas, Esgas 

We define Esgas as the energy necessary to carry the 
gas of a given galactic region outtside the galaxy. 

With reference to a spheric shell region of tickness 
ARi and internal and external radii Ri and Ri+i = 
Ri + ARi respectively, the energy Esgas can be written 
as follow: 

rRi+ARi 

E B9 as = I dL(r) (7) 



-s 



where dL(r) is the energy required to carry out a quan- 
tity dm = 47rr 2 f g (r, t)dr of gas. Here / 9 (r, i) is the den- 
sity function of gas, depending on the distance r from 
the galactic center and on the galactic age and t. There- 
fore, we have for dL(r): 



dL(r) 



dr'f(r') 



(8) 



where f(r') is the force between dm and the total mass 
(dark and luminous) within r' which reads as: 

GM(r')dm 



The range of integration in (8) can be reasonably lim- 
ited to an upper effective radius, since we only require 
that the gas leaves the galaxy. For this reason we re- 
place oo with (1 + a) R e , where R e is the effective 
galaxy radius and a > 0. As a consequence of this, 
dL(r) — J^ 1+a ^ Re dr'f(r'). In practice, we have ana- 
lyzed two cases: i) a = 0.5 which corresponds to the 
case in which the gas is carried out at 1.5i? e , and ii) 
a — oo which corresponds to the case in which the gas 
is carried out to infinity. The influence of a finite value 
of a will be discussed in section 3. 

The total mass at radius r', M(r'), is given by: 



M(r') 



Mr 



■ c (r') + M Lum (r') 



(10) 



where the luminous mass can be written as the sum of 
the gaseous and stellar components: 



M Lum (r') = M gas {r') + M.(r') 



(11) 



Assuming the same distribution for the two components, 
we have: 



M gas (r') = p g Fi(r'); 
M,(r') = p.F,(r'); 



f'T 

F l {r') = 4-n: 

Jo 



dsfi(s)s ; 



(12) 
(13) 

(14) 



The quantities appearing in equations (12), (13) and 
(14) can be determined on the basis of prescriptions 
available in the literature. Namely, we have for F;(r) 
(Jaffe 1983) 

r/r. 



Fi(r) tx 



(15) 



fir') 



(9) 



1 + r/r. 

with r. = 7i e /0.763, and for the dark matter distribution 
we have (Bertin Saglia and Stiavelli, 1992) 

M Dark (r) = 2 -Mf± { arct 9 {r/r ) 1 (16) 

with ro = 0.45i?_D ar fc. In the following, three different 
cases for the distribution of dark matter will be ana- 
lyzed: 

(i) M Dark = 

(ii) M D ark = WM L um and RDark = 10-R e 

(iii) M Dark = WM Lum and R Dark = R e 

These cases have been chosen since they represent ex- 
treme cases going from no dark matter (except the stel- 
lar remnants) to an amount of dark matter ten times 
the luminous one either distributed in a diffuse halo or 
tracing the luminous material. All the other possibilities 
will therefore be intermediate between these cases. 

In the case a = oo, Ri — 0, Ri+i = R e (i.e. one- 
zone model), considering luminous matter uniformly 
distributed and in absence of dark matter, eq. (7) gives: 
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3 G 

E Bga s = -^-M gas {2M f + M 3as ) (17) 

2.6 The thermal energy of the gas due to 
supernovae, E thsN 

The total thermal energy in the gas at the time t is: 

Eth SN {t)= [ €th SN {t-t')R S N{t')dt' (18) 
Jo 

where RsN{t) is the SN rate and eth SN (^sjv) is the ther- 
mal energy content inside a supernova remnant of both 
types, with tsN = t — t . e t h SN (tsN) evolves, according 
to Cox (1972), in the following way: 

e t h SN (tsN) = 7.210 50 e ergs (19) 

for < tsN < t c , and 

e thsN (t S N) = 2.2W 50 e {tsN/t c )- 0S2 ergs (20) 
for tsN > t c 

where e is the initial blast wave energy in units 
of 10 51 ergs, tsN is the time elapsed since the super- 
nova explosion and t c is the cooling time of a supernova 
remnant, namely 

t c = 5.710 4 e * /17 n - 9/17 yr (21) 

and no(t) = p(t)/ < m > is the average number den- 
sity of the ISM. Gibson (1997) adopted a more sophisti- 
cated treatment for the calculation of E t h SN , based on 
calculations by Cioffi McKee and Bertschinger (1988). 
However, given the uncertainties present in this kind of 
calculations we have decided to adopt the old formu- 
lation, common to most of the previous papers on the 
chemical evolution of ellipticals. We did not consider 
the energy injected into the ISM by stellar winds, since 
for massive galaxies is negligible, as shown by Gibson 
(1994). 

3 RESULTS AND CONCLUSIONS 
3.1 Model parameters 

We have calculated several chemical evolution models 
by varying the amount and distribution of dark matter 
as well as the parameter a and we have also considered 
either the case of a constant IMF in all the regions of 
the galaxy or the case of an IMF varying with the galac- 
tocentric distance. In Table 1 we show the choice of the 
parameters for all the models referring to the case with 
constant IMF. In particular, in this table we indicate the 
amount of dark matter (column two) and the distribu- 
tion of dark matter (column three). Finally, in column 
four we show the adopted value of the parameter a. We 



have used two different values for the IMF slope in the 
models with an IMF costant in radius: x = 0.95 and 
x = 1.35 (Salpeter 1955). 

In the framework of the models with a variable IMF 
we considered only one case in which all the parameters 
are the same as in the B2 model, except for the IMF 
slope. We will refer to this last case as the model IV, 
where the IMF slope varies from x = 0.95 in the in- 
ternal regions to x = 1.21 in the external ones. This 
choice is made in order to have the [Mg/Fe] constant 
in radius as it will be shown. Observational results seem, 
in fact, to indicate that [Mg/Fe] as a function of the 
galactocentric distance is roughly constant, at variance 
with what happens for the central [Mg/Fe] as a func- 
tion of galactic mass (Worthey et al. 1992; Carollo et 
al. 1993; Davies et al.1993). However, there is a warn- 
ing here, namely the trends on real abundances inferred 
from metallicity indices may be deceiving as recently 
suggested by Tantalo Bressan and Chiosi, (1997) and 
Matteucci and Ponzone (1997). In all the cases the val- 
ues of taw, the iron abundances and the [Mg/Fe] in 
the dominant stellar population will be calculated for a 
galaxy of 10 11 Mq of luminous mass and R e = 3Kpc. 

3.2 Galactic winds 

The results concerning galactic winds are summarized 
in Fig. 1 and Tables 2-14. Fig. 1 shows the behaviour of 
EBgas (solid line) and E t h SN (dashed line) for the 6 one- 
zone models. The occurence of galactic winds is affected 
by the presence of dark matter depending on the value 
of a. The value of a sets the distance at which the gas is 
ejected. Small values of a imply the presence of gas near 
the galaxy whereas large values of a imply that the gas 
is carried outside the gravitational field of the galaxy. 
This formulation allows us to study different cases unlike 
the model of Matteucci (1992) which corresponded only 
to the case of a small a. Although a dynamical model 
would be necessary to assess if the gas mass is definitely 
lost from the galaxy or if it remains around it, different 
values of the parameter a may influence the chemical 
enrichment of the ICM, in the sense that large values 
of a will predict a higher contribution from the galaxies 
to the gas of the ICM. To be more precise the influence 
of dark matter is very strong for all values of a when it 
is distributed as the luminous one (models CI and C2). 
However, when dark matter is distributed in a diffuse 
halo, its effect is relevant only for high values of a (Bl): 
in fact, in this case, we are considering the potential 
energy of the gas as the energy to carry the gas out the 
halo. On the other hand, a diffuse halo of dark matter 
does not affect the occurence of galactic wind if we use 
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small values of a (B2). This last result is in agreement 
with Matteucci (1992): she defined the potential energy 
as the energy to carry the gas out the galaxy, namely 
at the distance R e (a = 0). 

Tables 2-14 give the values of tGW for the models in 
the case of 10 zones with IMF slope equal to 0.95 (tab. 
2-7), 1.35 (tab. 8-13) and variable (tab. 14) (the zone 
tickness are obtained directly from data observation by 
Davies et al.,1993). 

/.From these tables we note that galactic winds are 
present only in a small number of outer regions for mod- 
els C. Namely, the wind occurs only for r > 1.5Kpc and 
r > l.OKpc for models CI and C2 respectively. Conse- 
quently, models C predict much smaller gas masses to 
be ejected from the galaxy and masses to form the hot 
coronae smaller than the those estimated from the X-ray 
emission (Forman Jones and Tucker, 1985; Canizares 
Fabbiano and Trinchieri, 1987). Such galaxies should 
also show active star formation at intermediate redshifts 
(z < 1), which seems not to be the case for normal el- 
lipticals (Matteucci 1994; Renzini 1995; Ziegler 1997). 
Therefore, we are tempted to reject the hypothesis of 
dark matter distributed like the luminous one (as al- 
ready pointed out by Matteucci 1992). 

With regards to models Al, A2, B2, we note they 
are very similar to those in the case of one-zone. For 
these models, galactic winds occur in the internal re- 
gions after several Gyr since galaxy formation (note that 
in the most internal region, with r < 30pc, the wind 
never occurs), while in the most external regions they 
occur before 1.0 Gyr from the epoch of galaxy forma- 
tion. 

This time spread of the galactic wind occurence 
with the galactocentric distance, has very important 
consequences on the chemical evolution in clusters of 
galaxies. Previous one zone models, in fact, including 
massive, diffuse, dark matter component (B2 models), 
undergo only a global "impulsive" ejection of their ISMs 
at tGW (the same result was pointed out by MG), also 
for different prescriptions of the IMF and SFR). There- 
fore, these one-zone models do not predict any chemical 
evolution in the ICM after nearly lGyr from galaxy for- 
mation, when galactic winds have already occurred in all 
the galaxies. On the other hand, the multi-zone mod- 
els predict chemical evolution in the ICM also at the 
present time because, even if there is only one global 
"impulsive" ejection of gas in all the regions of the 
galaxy, the most internal regions eject the processed gas 
very late. Another consequence, which we will see in the 
next paragraphs is the predicted behaviour of [Mg/Fe] 
in stars which tend to increase with increasing radius, 



unless a variable IMF is assumed. This could be a po- 
tential problem since the galactic wind occurring late in 
the most internal regions of ellipticals might mean active 
star formation at small red-shifts, which is not observed 
in these galaxies. However, as it is shown in Fig. 2 the 
SN II rates (and therefore the SFR) predicted by our 
best model IV starts to be negligible after t > 3 Gyr. 
This is due to the fact that the gas is consumed by star 
formation much before the galactic wind occurs. We do 
not show the SN rates for our model B2(x=0.95) since 
they are very similar to the results of model IV. 

3.3 Metallicity in stars 

According to the model of chemical evolution used, each 
region of elliptical galaxies has a composite stellar popu- 
lation containing stars of different ages and metallicities. 
This is, of course, true also in reality since the process 
of star formation is not instantaneous. In order to com- 
pare theoretical results with the metallicity abundances 
indicated by the observed indices we should first cal- 
culate the average iron and magnesium abundances in 
the composite stellar populations. The mass-averaged 
metallicity of a composite stellar population is defined 
following Pagel and Patchett (1975) as: 

1 f Sl 

< Z > M = tt / Z{S)dS (22) 
bl Jo 

where the subscript 1 indicates a specific time ti, which 
in our case is the present time, assumed to be 15 Gyr 
and S is the total mass of stars ever born. Accord- 
ing to the definition above we can calculate the aver- 
age iron and magnesium abundances, since the model 
provides the functions Xfs(S) and Xmq(S). However, 
it should be noted that, although the mass-averaged 
metallicity represents the real average metallicity of the 
composite stellar population, what should be compared 
with the metallicities deduced from observations is the 
luminosity-averaged metallicity. In fact, the metallicity 
indices, measured from integrated spectra, represent the 
metallicity of the stellar population which predominates 
in the visual light. The transformation from metallicity 
indices to real abundances should then be performed 
through either theoretical or empirical calibrations re- 
lating metallicity indices to real abundances. Galactic 
chemical evolution models (Arimoto and Yoshii, 1987; 
Matteucci and Tornambe, 1987) have shown that the 
two averaged metallicities are different, since metal poor 
giants tend to predominate in the visual luminosity so 
that, in general, the luminosity- averaged metallicity is 
smaller than the mass-averaged one. However, the dif- 
ference is negligible for galaxies with mass M > 10 11 M Q 
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which is the case we discuss here, and we are inter- 
ested mainly in studying metallicity gradients which are 
clearly not affected by the absolute value of the metal- 
licity. 

The results concerning the abundance of Fe in the 
dominant stellar population and the indices are sum- 
marized in the 3-d and 4-st columns of tab. 2-14. In 
these same tables we report the logarithm of (■§-) (col- 
umn one) and the time for the onset of galactic winds 
expressed in Gyr (column two) . In all the tables but ta- 
ble 14 column 5 shows the [^pf] ratio in the dominant 
stellar population. In table 14 instead we show in col- 
umn five the adopted values for the IMF slope as func- 
tions of the galactocentric distance. In this model [jQ] 
is constant as a function of the galactocentric distance 
([Mg/Fe] — 0.15), since we selected the IMF slope in 
order to achieve this result. The [Fe/H] of the dominant 
stellar population (column 3) is transformed into Mg2 
(column 4) by using the calibration of Worthey (1994) 
corresponding to an age of the dominant stellar popu- 
lation of 17 Gyr and [Mg/Fe]=0. 

In models A2 and B2 with x = 1.35 we obtain for 
[^] values of 0.55 and 0.025 respectively in the central 
and external regions. Correspondingly we obtain a gra- 
dient of A ^g^ = 0.27. We tried to use also other cal- 
ibrations available in the literature (Buzzoni Gariboldi 
and Mantegazza, 1992; Mould 1978; Barbuy, 1994) and 
we found differences either in the absolute values of Mg% 
up to ~ 0.1 dex. 

In fig. 3 the predicted values of Mg2 versus radial 
distance from galactic center are shown together with 
data from Davies et al. (1993). The data refer to galax- 
ies of the same size as our typical galaxy. In the figure 
the solid line represents the predictions of model IV, 
whereas dotted and dashed lines represent the predic- 
tions of model B2 with x = 0.95 and x = 1.35, respec- 
tively. Theoretical values are sistematically larger than 
observational data. This fact can be due to the uncer- 
tainties related to the metallicity calibration and to the 
stellar yields. On the basis of the previous remark it is 
interesting to note that a better matching between the- 
oretical and experimental data can be simply obtained 
by shifting the theoretical curves by a suitable amount. 
The shifts we applied are -0.079 dex for model IV, - 
0.084 for model B2 (0.95) and -0.05 for model B2(1.35) 
These shifts are well inside the differences we obtained 
by changing calibrations and/or by using models of dif- 
ferent ages as it is shown in fig. 4 where the predictions 
of our best model (model IV) are transformed according 
to the calibratrtions of Worthey (1994) for different ages 
of the dominant stellar populations and [Mg/Fe] =0. In 



any case, we do not think that is very important to 
fit the absolute values of Mg2 giving the still existing 
uncertainties, but that is more meaningful to fit the gra- 
dient itself and our models produce a very good fit. In 
particular, the model with x = 1.35 shows a metallic- 
ity gradient in very good agreement with data , but the 
value of [Mg/Fe] in the central region is —0.3, namely 
too low. In fact, the values suggested by population syn- 
thesis studies (Worthey et al. 1992; Weiss, Peletier and 
Matteucci, 1995) indicate positive values for this ratio. 
This result is mostly due to the late occurrence of galac- 
tic winds in the center of the galaxy which allow the 
supernovae of type la to restore the bulk of iron before 
star formation stops. However, as discussed by Gibson 
(1997), yields different than those adopted here may in- 
crease the value of [Mg/Fe] but only on short timescales, 
namely when the contribution of massive stars is still 
dominant. In fact at late times, when the SNe la are 
appearing, the [Mg/Fe] ratio decreases because of the 
newly injected iron. In addition, this model predicts that 
[Mg/Fe] should increase with increasing radius, clearly 
at variance with what inferred from metallicity indices. 
On the other hand, model B2 with x = 0.95 shows a 
metallicity gradient generally flatter than the observa- 
tional data but a better ratio [Mg/Fe] (+0.15 in the 
central region). However, the best model seems to be 
the model IV with a variable IMF which shows a good 
metallicity gradient and a correct value of [Mg/Fe] in 
the center and in the other galactic regions. This best 
model predicts also SN rates (expressed in units of SNu, 
namely SNe century^ 1 10~ 10 Lb ) and blue luminosity 
at the present time in good agreement with observa- 
tions (Turatto, Cappellaro and Benetti, 1994), namely 

(RsNIa)centre = 1.9 SNu, (R.S N 1 1) centre = SNu and 

(R S Nia) Rc = 0.17 SNu, (Rsnii)b^ = SNu. The pre- 
dicted integrated blue luminosity for the whole galaxy, 
calculated as in Matteucci (1994), is L© = 7 • 1O 9 Lb . 



4 CONCLUSIONS AND DISCUSSION 

In this paper we have explored the possibility of explain- 
ing metallicity gradients in giant ellipticals by means of 
models of chemical evolution including galactic winds. 
We have analyzed the case of an elliptical galaxy with 
initial luminous mass Mi um = lO n M , Re = 3 Kpc 
and with a dark matter halo ten times more massive 
than the luminous mass. 

We found that, in order to reproduce most of the 
observational constraints relative to elliptical galaxies: 

- the dark matter should not be distributed like 
the luminous matter but it should be more diffuse. This 
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result reinforces a previous analysis made by Matteucci 
(1992). 

-Abundance gradients in the stellar populations of 
ellipticals can be obtained by assuming that galactic 
winds develop first in the outer and later in the inner 
galactic regions, in agreement with the observational ev- 
idence of a correlation between the gas escape velocity 
and metallicity (Franx and Illingworth, 1990; Davies, 
Sadler and Peletier, 1993; Carollo and Danziger, 1994). 
The physical justification for this resides in the fact that 
the force required to carry out mass from the more ex- 
ternal regions is lower than that required to carry out 
mass in the more internal ones. This mechanism for the 
formation of gradients is different from the mechanism 
often invoked (Larson, 1974; Carlbergh 1984; Greggio 
1997). In this framework, an abundance gradient should 
arise because the stars form everywhere in a collapsing 
clouds and then remain in orbit with a little inward 
motion whereas the gas sinks further in because there 
is dissipation. This gas contains the metals ejected by 
evolving stars so that an abundance gradient develops 
in the gas. As stars continue to form their composition 
reflect the gaseous abundance gradient. 

Our model does not exclude a mechanism of this 
type which can be present together with the galactic 
winds occurring at different epochs in different regions. 
What we want to show here is that the differential occur- 
rence of galactic winds alone can explain the observed 
abundance gradients. On the other hand, a mechanism 
of inside-out galaxy formation, as that required to re- 
produce abundance gradients in the disk of our Galaxy 
(Matteucci and Francois, 1989: Chiappini, Matteucci 
and Gratton 1997) , does not seem to be suitable for 
elliptical galaxies. In fact, in order to reproduce the 
abundance gradient in the disk of our Galaxy one needs 
a large delay between the formation, due to infall of 
gas, of the more internal and the more external parts 
of the disk, of the order of several billion years (Chiap- 
pini et al., 1997). This would imply that most ellipticals 
should be still forming now with active star formation 
like in the disks of spirals, at variance with observations. 
The advantage of the model proposed here over other 
models such as that of Larson (1974) is that it explains 
the observed correlation between colors/metallicity and 
escape velocity 

A final remark concerns the star formation effi- 
ciency which we have assumed to be constant with 
galactocentric distance. We have also tried to vary this 
efficiency with galactocentric distance in the sense of 
having a higher efficiency in the central regions. We 
found that for moderate variations of the efficiency 



of star formation the abundance gradients steepen to- 
gether with gradient in [Mg/Fe]. On the other hand, if 
the variation of v is too strong the trend for the galactic 
wind is reversed, occurring later in the outermost re- 
gions with the consequence of cancelling the abundance 
gradients. 

- However, a potential problem with the mechanism 
for the formation of gradients suggested here is that it 
would predict an increasing [Mg/Fe] as a function of the 
galactocentric distance, at variance with observations, 
unless a variable IMF is assumed at the same time. 

- In order to obtain an almost constant [Mg/Fe] as 
a function of the galactocentric distance, as suggested 
by observations (Worthey et al., 1992; Carollo et al., 
1993; Carollo and Danziger 1994), one has also to as- 
sume that the IMF varies with radius. In particular, x 
should vary from 0.95 at the galaxy center up to 1.21 
in the outermost regions. This would imply a slightly 
higher percentage of massive stars in the innermost re- 
gions thus furtherly favoring the abundance gradients. 
However, the observed trend of [Mg/Fe] is inferred from 
the metallicity indices which are also affected by other 
parameters such as age. Therefore, it is very important 
to have reliable [Mg/Fe] ratios as a function of the galac- 
tocentric distances in ellipticals in order to assess this 
point. 

-The required IMF, in order to obtain positive val- 
ues of [Mg/Fe], especially in the galaxy center, as sug- 
gested by the observations (Worthey et al., 1992; Weiss 
et al., 1995), must have a slope x = 0.95 in agree- 
ment with previous studies (MG; Gibson and Matteucci, 
1997), although a variation in the stellar yields could 
achieve the same result (Gibson 1997) but only if the 
timescale for the galactic wind is short enough (Weiss, 
Peletier and Matteucci, 1995). 

-Another potential problem could be the fact that 
in our best model (model IV) the galactic wind does 
not develop until 12 Gyr in the most internal regions 
(< 150 pc). However, it should be noted that due to the 
fast consumption of gas star formation occurs at such a 
low level at 12 Gyr that would be hardly detectable. 

-Finally, we want to discuss the consequences of 
our model on the chemical evolution of the intracluster 
medium (ICM) although more details will be found in a 
future paper (Colafrancesco S., Matteucci F., Martinelli 
A., and Vittorio N., in preparation). The late occurrence 
of galactic winds in the more internal regions of ellip- 
ticals, if true, will have as a consequence a continuous 
chemical evolution of the ICM, at variance with the pre- 
dictions by models with only early winds (MG; Gibson 
and Matteucci, 1997) which predict no evolution of the 
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ICM starting from very high redshift (Z< 4 — 5). This 
prediction can be checked in the future with observa- 
tions of abundances in high red-shift clusters; prelimi- 
nary results from ASCA (Mushotsky and Lowenstein, 
1997) seem to indicate no evolution in the abundance 
of Fe for Z > 0.3. 
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B2 
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CI 
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C2 
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Table 1. Models considered. The parameter a, in 3-d col- 
umn, fixes the distance at which the gas is carried out the 
galaxy, namely (1 + a)R e 
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-1.385 
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0.18 
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0.59 


0.385 


0.20 


-0.424 


1.5 


0.55 


0.378 


0.22 


-0.279 


1.1 


0.52 


0.373 


0.24 


-0.137 


0.82 


0.48 


0.366 


0.25 


0.000 


0.62 


0.43 


0.357 


0.27 



Table 2. Model Al with the IMF slope (x) constant in radius 
and equal to 0.95 

FIGURE CAPTIONS 

Figure 1 Esgas (solid line) and E t h SN (dashed line) 
for the 6 models in the case of one-zone. 
Figure 2 The type II (dashed line) and la (continuous 
line) SN rates, expressed in units of SNe yr^ 1 , as a func- 
tion of time (expressed in Gyr) as predicted by model 
IV. The upper panel refers to the central region whereas 
the lower panel to the more external one (corresponding 
to an effective radius). 

Figure 3 Predicted and observed gradients of the 
metallicity index Mg 2 . The data refer to galaxies in the 
sample observed by Davies et al. (1993). The contin- 
uos line refers to model IV, the dashed line to model 
B2 (x=1.35) and the dotted line to model B2 (x=0.95). 
The error bars in the data are also shown. The calibra- 
tion adopted to transform [F/H] into Mg 2 is the one of 
Worthcy (1994) relative to an age of 17 Gyr. 
Figure 4 Gradients of the index Mg 2 as predicted 
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Table 4. Model Bl with x = 0.95 

by our model IV by using the calibrations of Worthey 
(1994) referring to different ages of the dominant stellar 
population. 
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Table 6. Model CI with x = 0.95 



>/ metallicity gradients in elliptical galaxies 

[hbt] 



[hbt] 



Log(^-) 


tow 




Mg 2 


[Mjl 
I Fe 1 


-2.233 


1 


0.68 


0.400 


0.15 


-1.929 


1 


0.68 


0.400 


0.15 


-1.385 


1 


0.68 


0.400 


0.15 


-0.963 


1 


0.68 


0.400 


0.15 


-0.824 


1 


0.68 


0.400 


0.15 


-0.602 


1 


0.68 


0.400 


0.15 


-0.424 


9.5 


0.67 


0.398 


0.15 


-0.279 


5.0 


0.65 


0.394 


0.17 


-0.137 


2.7 


0.61 


0.388 


0.19 


0.000 


1.4 


0.54 


0.376 


0.23 


* 7. Model C2 with x 


= 0.95 






Log(-g-) 


tGW 


m 


Mg 2 


[Mji 
L Fe ' 


-2.233 


1 


0.55 


0.378 


-0.32 


-1.929 


1 


0.55 


0.378 


-0.32 


-1.385 


1 


0.55 


0.378 


-0.32 


-0.963 


5.3 


0.50 


0.369 


-0.29 


-0.824 


3.6 


0.46 


0.362 


-0.27 


-0.602 


2.0 


0.35 


0.343 


-0.19 


-0.424 


1.26 


0.26 


0.329 


-0.13 


-0.279 


0.90 


0.20 


0.318 


-0.084 


-0.137 


0.72 


0.14 


0.308 


-0.044 


0.000 


0.56 


0.071 


0.296 


0.0030 



Table 8. Model Al with x = 1.35 
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Table 10. Model Bl with x = 1.35 
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0.071 


0.296 


0.0030 


0.000 


0.43 


0.025 


0.288 


0.030 


11. Model B2 with 


x = 1.35 






Log{£) 


taw 


m 


Mg 2 


[Mai 

1- Fe ' 


-2.233 


1 


0.55 


0.378 


-0.32 


-1.929 


1 


0.55 


0.378 


-0.32 


-1.385 


1 


0.55 


0.378 


-0.32 


-0.963 


1 


0.55 


0.378 


-0.32 


-0.824 


1 


0.55 


0.378 


-0.32 


-0.602 


1 


0.55 


0.378 


-0.32 


-0.424 


1 


0.55 


0.378 


-0.32 


-0.279 


11.9 


0.54 


0.377 


-0.32 


-0.137 


7.1 


0.52 


0.373 


-0.31 


0.000 


3.9 


0.47 


0.364 


-0.27 



Table 12. Model CI with x = 1.35 
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[hbt] 



Log(-Z-) 


taw 


m 


Mg 2 


L Fe > 


-2.233 


1 


0.55 


0.378 


-0.32 


-1.929 


1 


0.55 


0.378 


-0.32 


-1.385 


1 


0.55 


0.378 


-0.32 


-0.963 


1 


0.55 


0.378 


-0.32 


-0.824 


1 


0.55 


0.378 


-0.32 


-0.602 


1 


0.55 


0.378 


-0.32 


-0.424 


11.3 


0.54 


0.376 


-0.32 


-0.279 


5.2 


0.50 


0.369 


-0.29 


-0.137 


2.5 


0.40 


0.352 


-0.22 


0.000 


1.1 


0.25 


0.327 


-0.12 



Table 13. Model C2 with x = 1.35 



[hbt] 



Log(-g-) 


tGW 




Mg 2 


X 


-2.233 


1 


0.68 


0.400 


0.95 


-1.929 


1 


0.68 


0.400 


0.95 


-1.385 


12.3 


0.68 


0.400 


0.95 


-0.963 


4.4 


0.64 


0.393 


0.97 


-0.824 


3.1 


0.61 


0.389 


0.98 


-0.602 


1.8 


0.54 


0.376 


1.02 


-0.424 


1.2 


0.47 


0.364 


1.05 


-0.279 


0.80 


0.42 


0.350 


1.09 


-0.137 


0.58 


0.39 


0.333 


1.13 


0.000 


0.48 


0.17 


0.312 


1.21 



Table 14. Results related to 10-zones IV model. The IMF 
is radius dependent in order to have [Mg/Fe] constant (= 
0.15). 
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